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Abstract 

The  propagation  of  oscillatory  vaves  through  periodic  elastic  com- 
posites has  been  anailyzed  on  the  basis  of  Floquet  theory.  This  leads  to 
self-adjoint  differential  equation  systems  vhich  it  has  proved  conveni- 
ent to  solve  by  V8u:iational  methods.  Many  composites,  such  as  the  light- 
weight high-strength  boron-epoxy  material,  consist  of  strong  reinforcing 
components  in  a plastic  matrix.  The  latter  can  exhibit  iriscoelas^  ic  pro- 
perties which  can  have  a significant  influence  on  wave  propagation  charac- 
teristics. 

Replacement  of  the  elastic  constant  by  the  viscoelastic  complex 
modulus  changes  the  mathematical  structixre  so  that  the  differential  equa- 
tion system  is  no  longer  self-adjoint.  However,  a modification  of  the 
variationsil  principles  is  suggested  which  retains  formal  self-adjointness, 
and  yields  variational  principles  which  contain  additional  boundary  terms. 
These  are  applied  to  the  determination  of  wave  speeds  and  mode  shapes  for 
a laminated  composite  made  of  homogeneous  elastic  reinforcing  plates  in 
a homogeneous  viscoelastic  matrix  for  plane  waves  propagating  normally 
to  the  reinforcing  plates.  These  results  agree  well  with  the  exact  solu- 
tion which  can  be  evaluated  in  this  simple  case.  The  variational  principles 
permit  solutions  for  periodic,  but  otherwise  arbitrary  variation  of  ma- 
terial properties. 


■r 


INTRODUCTION 


4 » 


A composite  medium  typically  consists  of  a matrix  material  with  an 
embedded  reinforcing  material  in  the  form  of  fibers  or  laminations.  A 
laminated  ccxnposite  cein  usually  be  modelled  as  alternating  layers  of 
matrix  and  filament  arranged  in  a periodic  manner,  while  a fiber  rein- 
forced composite  can  often  be  represented  by  a homogenious  matrix  in 

1 

which  a two  dimensional  doubly  periodic  array  of  filaments  is  embedded.  , 

A convenient  method  for  analysing  wave  propagation  in  composite  j 

media  is  the  use  of  Ploquet  theory.  This  approach  has  been  recently 

i 

used  by  several  authors  [1-7]  to  study  steady  state  wave  propagation  in 
periodic  elastic  con5>osites.  Of  these,  some  [1-5 ] used  veuriational  methods 
while  others  [6-7]  used  direct  numerical  methods  involving  discretiza- 
tion of  the  governing  differential  equations  and  associated  quasi-peri- 
odic  boundary  conditions.  | 

In  practice,  the  matrix  in  a composite  is  often  a polymer  which  exhi- 
bits viscoelastic  properties.  This  leads  to  dissipation  as  well  as  dis- 
persion. Free  wave  propagation  in  an  infinite  laminated  viscoelastic 

canposite  was  studied  by  the  present  authors  in  [8]  using  Ploquet  theory  | 

and  finite  difference  methods.  In  the  viscoelastic  case,  forced  steady 
wave  propagation  leads  to  real  frequencies  and  complex  wave  numbers 
(attenuating  waves)  while  free  waves  with  no  applied  tractions  have  com- 
plex frequenoies  and  real  wave  numbers  (damped  waves)  [8]. 

Consideration  of  viscoelastic  properties  of  the  matrix  leads  to 
complex  viscoelastic  moduli  which  are  functions  of  the  frequency  that 
can  be  real  or  complex.  EJven  in  the  elastic  case  [l]  iu  is  convenient 
to  utilize  complex  analysis  to  incorporate  the  changing  phase  as  the 
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wave  traverses  the  periodic  medium,  so  that  the  displacement  is  expressed 

by 

- / V i(ut 

u = u(x)e 

where  u is  the  frequency,  t is  time  and  u(x)  is  complex.  Applica- 
tion of  variational  methods  involves  integrals  of  the  type 

, * 

a/2  du  du  ^ 

Kuj.Uj)  = p( x)« 

where  ri(x)  is  the  elastic  moduJ.us  veiriation,  p (x)  that  of  the  deisity, 
a is  the  periodic  cell  length  and  * denotes  the  complex  conjugate 
For  elo  ..icity  with  n e.nd  frequency  u real,  the  integrals  are 
Hermitian:  I(Uj^,U2)  = I (u2,Uj^)  and  self-adjoint  differential  operator 
systems  result  for  determining  the  displacement  or  stress  fields.  For 
viscoelasticity,  when  i and  o)  are  complex,  the  integrals  are  not 
Hermitian  since  n h and  the  resulting  differential  operators  are 
non-self-adjoint.  However,  if  one  formally  bases  the  theory  on  real 
analysis,  so  that  the  complex  conjugate  sign  is  removed  from  the  integral 
even  though  u(x)  is  complex,  the  integral  is  then  symmetric:  l(u^,U2)  = 
l(u2,Uj^),  and  hence  formally  self-adjoint  differential  operators  are  gen- 
erated. This  is  so,  with  u(x)  complex,  even  though  the  integral  l(u,u)  is 
then  no  longer  physically  an  enorgy  integral  and  the  variational  princi- 
ple no  longer  expresses  Hamilton's  principle.  This  formulation  introduces 
an  additional  boundsry  term,  so  that  the  system  is  formally  self  adjoint: 
a self-adjoint  (.nerator  but  not  with  th(>  appropriate  boundary  conditions. 
These  boundary  terms  can  be  included  to  generate  variational  principles, 
and  the  ones  presetted  here  are  formulated  in  this  way. 


- 3 - 


In  this  pap'-r,  the  strain  energy  [l],  compliinent8u:y  energy  [2]  and 
Hellinger-Reissner  [3-5l  variational  principles  sure  extended  so  that 
they  can  be  used  in  the  viscoelastic  case.  As  an  illustration,  the  ex- 
tended version  of  the  Hellinger-Reissner  [U]  vsuriational  principle  is 
applied  to  the  problem  of  free  waves  in  an  infinite  one  dimensional  visco- 
elastic composite  [8].  This  principle  ia  chosen  because  it  gave  very 
accurate  results  in  the  elastic  case  [U].  A composite  with  two  homogene- 
ous layers  per  cell  is  studied  with  the  filament  elastic  and  the  matrix 
modelled  as  a three  element  solid.  The  Rayleigh-Ritz  method  is  used  to 
obtain  dispersion  relations  and  mode  shapes  and  the  results  compared  with 
the  exact  solution  which  exists  in  this  case.  For  two  homogeneous  com- 
ponents the  elastic  solution  has  been  obtained  in  closed  form  [9»2J  and 
this  can  be  adapted  to  the  case  of  linear  viscoelastic  components.  The 
question  of  rapidity  of  convergence  to  the  exact  solution  is  discussed. 

The  present  method  can  be  used  to  study  wave  motions  in  general  one- 
dimensional periodic  viscoelastic  composites  which  otherwise  exhibit  arbitr- 
ary veiriations  of  material  properties.  Waves  in  fiber  reinforced  visco- 
elastic composites  can  be  studied  by  using  similar  variational  principle" 


in  two  dimensions. 


GOVERNING  EQUATIONS 


We  consider  wave  propagation  in  a one  dimensional  laminated  visco- 
elastic composite,  and,  in  peo'ticular,  choose  a two  material  composite 
as  shown  in  Fig.  1.  This  example  is  chosen  because  the  variational 
principles  to  be  presented  in  the  next  section  can  then  be  easily  com- 
P80*ed  with  those  in  the  elastic  case  [l-5]*  Variational  principles 
for  any  general  one-dimensional  inhomogenious  periodic  viscoelastic 
medium  can  be  easily  obtained  by  suitable  modification  of  these  princi- 
ples. 

The  composite  covers  the  f\ill  space  - oo  < x < <».  We  study  one- 
dimensional  strain  waves  propagating  in  a direction  x normal  to  the 
interface  planes.  For  harmonic  waves,  the  stress  and  displacement  are 
of  the  form  5 = o(x)e^^^  and  Q = u(x)e^*^^.  The  quantities  o and 
u are,  in  general,  complex. 

Within  each  cell,  the  governing  differential  equations  are 

^ + p (x)(o^u  = 0 (l) 

dx 

a = n(x,(u)^  (2) 

dx 

where  n = A + 2y  in  terms  of  complex  LamS  viscoelastic  moduli  of 
the  constituent  materials.  The  quantities  n and  p are  periodic  with 
period  a 

n(x+a;(i))  = Ti(x;(o)  , p(x+a)  = p (x) 

and  within  each  cell  they  arc:  discontinuous  functions  defined  by 
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R 


-a/2  < 

X < -b/2 

n2(“)»  Pg 

-b/2  < 

X < b/2 

n^(o)),  p^ 

b/2  < 

X < a/2 

E 


We  note  that  and  rtg  are,  in  general,  congjlex  functions  of  the 
possibly  ccxnplex  frequency  ui. 

By  Floquet  theory  [l]  the  displacement  and  stress  satisfy  the  quasi- 
periodic  boundary  conditions  across  each  cell 


i 

i 

f 

u(a/2) 

= u(-a/2)e^‘l^ 

(3) 

h 

1 

1 

o(a/2) 

= o(-a/2)e^*^^ 

(h) 

t 

f 

t 

where  q is  the 

wave  number. 

i 

The  interfaces  eire  assumed  perfectly  bonded  so  that  u and  n 

are 

s 

:■ 

'f 

i 

continuous  across 

them 

f 

i 

1 

“<v> 

II 

o 

X 

1 

o 

II 

± b/2 

(5) 

; 

r 

»<v> 

= a(x  -)  X = 

O 0 

± b/2  . 

(6) 

VARIATIONAL  PRINCIPLES 


1.  Strain  energy  principle 

Let  us  choose  a pair  of  sectionally  continuous  and  differentiable 
functions  u and  o which  satisfy  the  constitutive  equation  (2),  the 
quasi-periodic  boundary  conditions  (3)  and  (U)  and  continuity  condition 
(5).  The  solution  (u.o)  which  also  satisfies  the  equation  of  motion  (l) 
and  the  continuity  condition  (6)  is  given  by  the  variational  equation 

6 l(u)  + a(a/2)S  u(a/2) (l-e  = 0 (T) 

where 

t 2 2 2*^ 

Ku)  = / (-  ■ jd>  . 

-a/2 


Proof : 

Taking  variations  with  respect  to  u and  after  carrying  out  the 
necessary  integrations,  we  obtain 

6 I(u)  + o(a/2)6  u(a/2)(l-e“^^^®')  = J l^(n(x)|^)+p  (x)o)^u>5  udx 

-a/2^  ^ 


{cl 


du/ 


+ (o(a/2)-n^;^(a/2))6  u(a/2) 


+ n,  J^(-a/2)6u(-a/2)-o(a/2)5u{a/2)e' 

i OX 


+ < n(b/2)f^(b/2)6u(b/2)  > 
QX 


+ < n(-b/2)fi(-b/2)6u(-b/2) 
ax 


Here  we  use  the  notation 


< g(x  ) > = g(x^)  - g(x~) 

O O o 


- 6 


2iqa 
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If  u and  0 are  restricted  to  tne  subclass  of  functions  which 
satisfy  equations(2),  (3),  (^)  and  (5).  this  expression  reduces  to 

a/2 

J (•^  + p (x)w^u)  Ax  dx  ♦ < o(b/2)  > 6u(b/2)  ♦ < o(-b/2)  > 6u(-b/2)  . 

-a/2 

This  expression  vanishes  for  arbitrary  6xx  if  and  only  if  u and 
0 satisfy  equations  (l)  and  (6). 

2.  Coiaplenentary  energy  principle 

Let  us  choose  a pair  of  sectionally  continuous  and  differentiable 
functions  u and  o which  satisfy  the  eqixation  of  motion  (l),  the 
quasi-periodic  boundary  conditions  (3)  and  {h)  and  continuity  condition 
(6).  The  solution  (u^j ) which  also  satisfies  the  constitutive  equation 
(2)  and  the  continuity  condition  (5)  is  given  by  the  variational  equation 

6j(o)  + u(a/2)5o.'a/2){l-e‘^^^*)  * 0 (8) 

where 


and  the  prime  denotes  differentiation  with  respect  to  x. 

Proof: 

Taking  variations  with  respect  to  o and  after  carrying  out  the 
necessary  integrations  we  obtain 


a/2 

SJ(o)  » u(a/2)«o(a/2)[l-e-=^'>^l  = f - 


1 6a  dx 


+ lu(a/2)  + ■■  M 6a (a/2) 


- 6a(-a/2)-u(a/2)6a(a/2)e~^''*®^ 


p (b/2)o)‘' 


6a(b/2)  > 


g'(-b/2] 


p (-b/2)oj'^ 


6a(-b/2)  > . 


If  u euid  0 are  restricted  to  the  subclass  of  functions  which 
satisfy  equations  (l),  (3)»  (^)  and  (6),  this  expression  reduces  to 


/ ^ ^ u(b/2)>  5 o(W2)  + < u(-b/2)  > 5o(-b/2). 

-a/2 

This  expression  vanishes  for  arbitrary  6a  if  and  only  if  u and 
a satisfy  equations  (2)  and  (5). 

3.  Hellinger-Reisiner  principle. 

Let  us  choose  a pair  of  sectionally  continuous  and  differentiable 
functions  u and  a which  satisfy  the  quasi -periodic  boundary  condi- 
tions (3)  and  (^4)  and  continuity  condition  (5).  The  solution  (u,o)which 
also  satisfies  the  equation  of  motion  (l),  the  constitutive  equation  (2) 
and  continuity  condition  (6)  is  given  by  the  variational  equation 


6K(u,o)  + o(a/2)6u(a/2) (l-e  = o 
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a/2  2 2 2 

>aere  K(u,o)  = /("’§*  2^  * )■>» 

-a/2 


Proof 

Taking  variations  with  respect  to  u smd  o and  after  carrying 

out  the  necessary  integrations  we  obtain 

a/2  a/2 

6K(u,o)  +o(a/2)6u(a/2)(l-e"^^‘^^)  = ^ -fp  (x)u^u)6u  dx  + J 

-a/2  -a/2 


^)6o  dx 
dx 


+ o(-a/2) 6u(-a/2)  - o(a/2)6u(a/2)e 

+ < o(b/2)6u(b/2)  > + < o(-b/2)6u(-b/2)  > . 

If  u and  0 are  restricted  to  the  subclass  of  functions  which  satisfy 
equations  (3),  (^)  and  (5)»  this  expression  reduces  to 

a/2 


a/2 

/ 

-a/2 


du 


/ dx  t / (^-5; 


-a/2 


) 60  dx  + < o(b/2)  > 6u(b/2) 

+ < o(-b/2)  > 6u(-b/2). 


This  expression  vanishes  for  arbitrary  6u  and  60  if  and  only  if 
u and  0 satisfy  equations  (l),  (2)  and  (6). 

I4.  A general  variational  principle 
^he  variational  equation 

6[K(u,a)+L^+L  + o(a/2){ 6u(a/2)-6u(-a/2)e 


+ {u(a/2)e"^'^'' 


-u(-a/2))  (So (a/2)  = 0 


(10) 


where 


' \ - 
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= - |{o(b"/2)+o(b‘^/2)l  < u(b/2)  > 

L_^  = - ^a(-b"/2)+a(-bV2)]  < u(-b/2)  > 

and  K(u,o)  is  defined  in  equation  (9)  is  completely  equivalent  to  the 
original  boiindary-value  problem  (l)~{6). 

Proof. 

Taking  variations  with  respect  to  u and  o and  after  carrying 

out  the  necessary  integrations,  we  obtain 
a/2  a/2 

JT  (g  + p(x)A)6udx+  JT  (^-|^)«adx 
-a/2  -a/2 

+ {u(a/2)e  ^'^^-u(-a/2)}5o{a/2)  - {o(a/2)^'^^-a(-a/2)}  6u(-a/2) 

- i < u(b/2)  > {6a(bV2)+6o(b”/2)}  - ^ < u(-b/2)  > {6a(-bV2)+60(-b"/2)} 

+ I <o(b/2)  > {5u(b'*’/2)  + 6u(b'/2)}+  | < o(-b/2)  > {6u(-b'*‘/2)+6u(-b"/2)} 

This  expression  vanishes  for  arbitrary  6u  and  6o  if  and  only  if 
u and  0 satisfy  the  complete  set  of  equations  (l)-(6). 


A NUMERICAL  EXAMPLE 


1.  The  problem  an^  method  of  solution. 

As  an  illustration,  we  ccnsider  an  application  of  the  extended 
version  of  the  Hellinger-Reissner  principle  (9).  This  principle  offers 
greater  flexibility  than  the  strain  energy  or  complementazy  energy 
principles  and  Nanat-Nasser  [1»]  used  its  elastic  smeLLog  to  obtain  very 
accurate  results  for  wave  propagation  in  elastic  composites.  The  gen- 
eral variational  principle  (lO)  is  not  usM  since,  in  the  one-dimen- 
sional case,  it  is  easy  to  choose  test  functions  for  u and  o which 
are  continuous  and  satisfy  the  quasi-periodic  bovmdsry  conditions  (3) 
and  (U).  The  two-dimensional  version  of  (lO)  should  prove  very  useful 
for  study  of  wave  propagation  in  fiber  reinforced  viscoelastic  composites. 

The  exact  dispersion  relation  for  this  problem  was  obtained  by  Sun, 
Achenbach  and  Herrmann  [9] 

2 

cosQa  = cos  )cos(— ) - sin(“^^~^^)sin(— ) (11) 

Cl  Cj  2p  Cl  c^ 

2 2 2 
where  c^  = , c^  = tig/pg  and  p = 

The  exact  equation  for  the  mode  shape  is  given  in  [2]. 

We  note  that  since  we  are  studying  free  waves  in  an  infinite 
viscoelastic  composite  is  a complex  function  of  the  complex 
frequency  u.  The  wave  number  q is  real. 

We  obtain  approximate  solutions  for  various  numbers  of  terms 
in  series  expansions  for  u and  o and  compare  these  with  the  exact 
solution.  This  tests  the  efficiency  of  the  variational  principle  (9). 
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Ve  define  dioeneionless  quantities 


X = 

x/a 

, u(x) 

= u(x)/a  , S(x)  = a(x)/ng 

Q = 

qa 

. X(u,s) 

* K(u,a)/nga  , 

a = 

2 

aPi 

2, 

“ /hg  » 

0 = Pg/Pj^  » K = hg/n^  , 

a * 

b/a 

In  terns  of  these,  equation  (9)  becomes 


6K(U,S)  + S(l/ii)5U{l/2){l-e”^^*^)  = 0 


where 


1/2 


-1/2 


(12) 


Me  e]q>and  U and  S in  series 


^i(Q+2iTj)x 


J=-n 


-n 


* 2 v.<^>  = 2 


D e 
k 


i(Q+2nk)x 


k=-r 


-n 


(13) 


It  is  clear  that  U and  S satisfy  (3),  (*♦),  (5)  and  (6). 
Substituting  (13)  into  (12)  and  after  performing  the  necessary  inte- 


grations, ve  obtain 
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K(U,S)  = 


2 <<=j<i.^>jVjVw!^> 

J»k*-n 


u 

S(l/2)6U(l/2)(l-e"^^‘^)  = ^ J®^^5 


where 


1/2 


J ,k=-n 


sinQt(6-l)sin{a(Q-Hr(.1+]s 

2(Q+ir(j+k)) 


if  Q -Tt(j+k) 


[•>•0(6-3 

2 


if  Q = "iiCj-^k) 


- sin{g(ft»Ti(.1-*-k)  }}+^ 


s inQ-s  in  { a { Q+tt  ( j +h 


2(Q+Tr(j+k)) 


if  Q / -TT(j->k) 


if  Q = ~n(j+k) 


£ 
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dU  (X) 

I 

-1/2 


(-l)J'*'^i(Q»2ir.l)8ina 

Q+ir(j+k) 


i(Q+2irJ) 


(^_g-2iQj^i(Q+it(j+k)) 


if  Q 5^  -ir(j+k) 


if  Q = -ir(j+k) 


Taking  veiriations  and  equating  the  coefficients  of  fiC  and  6D 

d 

to  zero,  we  obtain  the  matrix  eqiiation 


u s 

(where  m = 2n+l).  The  elements  of  I and  I are  functions  of  the 

complex  frequency  ui  through  fi  and  n.  Thus,  for  ar»y  given  wave  number 
Q,  the  frequency  u can  be  obtained  by  solving  the  equation 


(15) 


det[N(o),q)]  = 0 


(16) 


where  N is  a 2m  x 2m  matrix 


» 1 

,su  , 

21  ' 

J -I 

— _ _ 

otS 

21 

Once  an  eigenvalue  u)  is  known,  the  correpsonding  eigenvector  and 
hence  the  mode  shapes  U and  S can  be  easily  calculated. 


(19) 
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We  now  rewrite  the  matrix  N(s)  in  order  to  show  the  frequency 
dependence  explicitly.  Thus,  equation  (l6)  becomes 


det 


As 


.2 


I 

u 


B 


-I 


I (1+s't  ) 

1 

I -L  t 


= 0 


(20) 


where  A,  B,  E,  F,  G are  (m  x m)  matrices  given  by 


A.,  - . IM)'’^‘‘slnS>(6-l)sln(arsaUaml  if  « ^ ..(j.k) 

n (Q+if(j+lc)) 


(H-g(e-l)) 

G 

n 


if  Q = -itCj+k) 


„ ..SU  ,3U 

Jk  jk  Jk 


su 


®Jk  * ^kj 


r - J1  [(-l)’^'*'^sinQ-sin{a(Q-»Tr(j-»k))}] 
■Jk  Q+Ti(j+k) 

= n^(l-a) 


n ^ sin{g(Q+Tr(j-t-k))} 
Jk  Q+Ti{J+k) 


= g 


if  Q -ir(j+k) 
if  Q = -Tt(j+k) 
if  Q ^ -TT(j+k) 

if  Q = -n(j+k) 


The  elements  of  A,  B,  E,  F and  G are  known  once  the  wave  number 
emd  the  other  parameters  of  the  pioblem  are  specified.  Equation  (20)  re- 
presents a nonlinear  eignevedue  problem  since  many  of  the  elements  of  N 
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4 

I 

i 

s 

a 

I; 

I 


are  nonllneso*  functions  of  s.  It  is  possible  to  reduce  this  to  a standard 
eigenvalue  problem  so  that  s*  is  an  eigenvalue  of  a new  matrix  M of  order 
3m  X 3m.  This  increases  the  size  of  the  problem  but  we  can  now  use  any 
of  the  standard  con5)uter  algorithms  for  determination  of  eigenvalues  of  a 
matrix. 

We  observe  that 


det 


»?11 


M, 


^12 


•■21  . “22 
I 


.-1 


if  the  matrices  are  conformable  and  exists,  so  that  equation  (20) 

.LX 

can  be  expanded  as 


where 


det[l  s'^-i-P  s’^-i-H  s'+yH]  = 0 


(21) 


It 

S“^(F-»G6,  ) 

[ 

H = 

V 

t 

1 

• 

Y = 

6 /t 
l'  1 

P 

R = 

t^(G•^F)  , 

and  I is  the  unit  matrix. 

Now  s*  is  an  eigenvalue  of  a 3m  x 3m  matrix  M defined  as 


M = 


0 ! I ! 0 

t- 

0 I 0 ' I 

- I ••  I ~ 



! -yH  • -h  1 -p 
L ' : '1  ' 


(22) 
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This  can  Te  easily  proved  by  expanding  the  equation 

det[*^I  s']  = 0 . (23) 

2.  Results  and  conclusions. 

Numerical  calculations  are  carried  out  for  the  following  values  of  the 
parameters  [8] 

G 

n = 4 and  50  0 = 3 

6l  = .70  t^  = .0^55  a = .5 

m = 7 and  11 

The  wave  number  Q is  specified  and  only  the  range  0 _<  Q ^ v needs 
to  be  studied  [lO]  because  the  Floquet  form  of  the  solution  is  not  uniquely 
determined. 

The  dispersion  ctirves  are  obtained  by  calculating  the  eigenvalues  of 
M (equation  (22))  using  an  IBM  360/67  computer.  A standard  QR  algorithm 
is  used.  Once  the  frequency  is  known  for  a certain  Q we  calculate  the 
corresponding  eigenvector  of  the  matrix  N (equations  15,  20).  Now  the 
mode  shape  U and  S are  obtained  from  equation  (13). 

Figures  2,  3,  and  5 show  the  real  and  imaginary  parts  of  the 
frequency  u as  a function  of  the  wave  number  Q for  the  first  5 modes. 
The  approximate  solutions  for  m = 7 and  11  are  compared  with  the  exact 
solution  obtained  from  (ll).  The  agreement  is  excellent  with  m = 11 
for  the  first  two  modes  but  the  approximate  solutions  become  progressively 
inaccurate  for  higher  modes.  The  maximum  error  in  the  real  part  of  the 
frequency  for  m = 11  is  of  the  order  of  10?J  and  occurs  for  the  fifth 
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mode  for  Q = it.  The  error  is  much  higher  for  the  imaiginary  ptirt  of  the 
frequency  for  the  fifth  mode.  We  recall  that  2,^  determines  the  rate 
of  exponentieul  decay  in  time. 

A comparison  with  [8]  shows  that  in  this  case  the  finite  difference 
scheme  is  more  efficient.  Also,  a comparison  with  [U]  shows  that  the 
Hellinger-Reissner  principle  does  not  work  as  efficiently  as  it  does  in 
the  elastic  case.  The  fundamental  mode,  however,  is  obtained  very  accur- 
ately with  m = 7 and  this  mode  dominates  since  the  higher  modes  decay 
faster . 

Mode  shapes  are  shown  in  Figs.  6 and  T«  The  real  parts  of  the  dis- 
placement ar*d  stress  are  shown  for  the  fundamental  mode  for  Q = ir/2. 

The  normalzing  constant  is  chosen  as  in  [8]  so  that  U(0)  = 1 + Oi.  The 
results  for  m = 11  are  in  excellent  agreement  with  the  exact  solution 
even  though  the  stress  and  displacement  are  expanded  in  series  of  smooth 
functions  (13)  while  the  solutions  contain  sharp  discontinuities  in  slope 
at  the  interfaces.  The  stress  solution  is  more  accurate  than  in  [8]  because 
here  it  is  obtained  independently  rather  than  being  obtained  from  U by 
numerical  differentiation. 

The  Hellinger-Reissner  principle  enables  us  to  obtain  continuity  of 
both  stress  and  displacement  across  the  interface  while  if  the  strain 
energy  principle  is  used  together  with  a smooth  displacement  field  a 
stress  discontinuity  results  [l]. 

To  s.um  up,  this  paper  presents  some  varational  principles  which  can 
be  used  to  study  wave  propagation  in  laminated  viscoelastic  composites.  A 
numerical  example  gives  encouraging  results.  These  principles  can  be 
easily  extended  to  higher  dimensions  to  study  wave  propagation  in  fiber 
reinforced  viscoelastic  composites. 


ACKNOWLEDGMENT 

This  research  was  supported  by  the  Office  of  Naval  Research  under 
Contract  NOOOIU75-C-O698  with  Stanford  University,  Steuiford,  California. 


- 20  - 


REFERENCES 


1.  Kohn,  J.A.  Krumhansl  and  E.H.  Lee,  Variational  methods  for  dis- 
persion relations  and  elastic  properties  of  composite  materials. 

J.  Appl.  Mech.  39,  327  (1972). 

2.  E.H.  Lee,  A survey  of  variational  methods  for  elastic  wave  propagation 
analysis  in  composites  with  periodic  structures,  I^amics  of  Composite 
Materials,  ASME  122  (1972). 

3.  S.  Hemat  Nasser,  General  veiriational  methods  for  waves  in  elastic 
composites,  Jn.  Elasticity  2,  73  (1972). 

1*.  S.  Nemat  Nasser,  Harmonic  waves  in  layered  composites.  J.  Appl. 

Mech.  39,  850  (1972). 

5.  S.  Nemat  Nasser  and  F.C.L.  Fu,  Harmonic  waves  in  layered  composites: 
Bounds  on  frequencies.  J.  Appl.  Mech.  Ul,  288  (197^) • 

6.  W.H.  Yang  and  E.H.  Lee,  Modal  analysis  of  Floquet  waves  in  composite 
material,  Jn,  Appl.  Mech.  Ul,  429  (1974). 

7.  G.H.  Golub  L.  Jenning  and  W.H.  Yang,  Waves  in  periodically  structured 
medium.  To  appear. 

8.  S.  Mukherjee  and  E.H.  Lee,  Dispersion  relations  and  mode  shapes  for 
waves  in  laminated  viscoelastic  composites  by  finite  difference 
methods.  Jn.  Computers  and  Structiu*es,  5,  279  (1975). 

9.  C.T.  Sun,  J.D.  Achenbach  and  G.  Herrmann,  Continuum  theory  for  a 
laminated  medium,  Jn.  Appl.  Mech.  35,  467  (1968). 

10.  E.H.  Lee  and  W.  Yang,  On  waves  in  composite  materials  with  periodic 
structure,  SIAM  Jn.  Appl.  Math.  25,  492  (1973). 


- 21  - 


I 


I 


I 


I 


I 


I 


Q 


EXACT  SOLUTION 

, I plane  waves 

7 PLANE  WAVES 


Figure  2.  Dispersion  Ciirves  u_  versus  Q 
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Figure  5.  Damping  Curves  versus  Q 
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Figure  6.  Displacement  and  Stress  Distribution  for  the  First  Mode 
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